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LONG-TERM  GOALS 

•  To  enhance  our  ability  to  detect  self-propelled  submerged  bodies  (SSBs),  e.g.  submarines  and 
AUVs  at  a  distance.  Little  work  has  been  done  to  improve  our  understanding  of  how  we  can 
sense  changes  in  the  environment  caused  by  internal  waves  (IWs)  and  wakes  generated  by  the 
passage  of  an  SSB. 

Simulations  of  self-propelled  submerged  bodies  have  been  conducted  in  idealized  environments  with 
simple  environmental  geometries  and  stratifications.  The  focus  of  those  simulations  was  on  improving 
the  submarine,  i.e.  reduced  drag,  less  vortex  production,  improved  maneuverability  and  greater  speed. 
Thus  many  of  the  simulations  focused  on  the  near  field  and  on  complex,  realistic  submarine  geometry. 
Some  simulations  looked  at  the  far  field  but  usually  for  simple  shapes  and  simple  background 
stratification. 

OBJECTIVES 

•  To  incorporate  the  ability  to  simulate  a  SSB  into  the  NRL-MIT  nonhydrostatic  model.  A  SSB 
can  be  simulated  as  a  smooth,  bluff  or  streamlined  body  moving  through  the  model  domain.  It 
can  also  be  simulated  as  a  point  momentum  source  with  or  without  swirl  (to  better  simulate  a 
propeller)  to  generate  the  extreme  far  field. 

These  simulations  will  focus  on  the  environmental  impact  of  a  SSB  to  improve  our  understanding  of 
the  signatures  that  would  be  measured  by  acoustical  or  E+M  systems.  To  this  end  we  are  interested  in 
hindcasting  realistic  coastal  ocean  environments  with  accurate  stratification,  currents,  ambient  IW 
fields  and  bathymetry  and  simulating  changes  in  this  environment  caused  by  the  wake  and  IWs 
generated  by  the  SSBs  movement.  In  a  stratified  fluid  the  wake  collapses  vertically  and  spreads 
horizontally.  Thus  its  signature  becomes  very  two  dimensional.  However,  the  collapse  also  generates 
IWs  which  propagate  three  dimensionally  thus  yielding  additional  potential  signatures  of  the  SSB. 
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APPROACH 


A  relatively  smooth  object  can  be  inserted  into  the  NRL-MIT  model  using  the  immersed  boundary 
method  (IMB).  Several  versions  of  this  technique  have  been  developed  e.g.  Tseng  and  Ferziger  (2003), 
Goldstein  et  al.  (1993),  Lai  and  Perskin  (2000),  and  Fadlun  et  al.  (2000).  Tseng  and  Ferziger 
demonstrated  their  version  by  using  it  to  simulate  a  complex  bottom  bathymetry  in  the  MIT  model. 

The  submerged  body  must  be  self-propelled  because  there  are  significant  differences  between  towed 
(drag)  wakes  and  self-propelled  (momentumless)  wakes.  The  details  of  the  geometries  and  relative 
motions  of  the  parts  of  the  SSB,  e.g.  sail,  diving  planes,  and  propeller(s)  or  jet,  are  critically  important 
in  determining  the  near  wake.  However,  much  of  the  detail  is  lost  in  the  far  wakes  particularly  after 
wake  collapse  at  approximately  7Nt  (where  N  is  the  Brunt-Vaisala  frequency).  However,  swirl  from 
the  propulsion  system  can  remain  a  factor  for  many  nondimensional  time  periods. 

We  plan  to  simulate  the  SSB  as  a  smooth,  self-propelled  body  with  and  without  swirl.  Since  we  are 
interested  in  the  far  wake  and  the  IW  fields  that  are  generated  from  the  wake,  particularly  during  wake 
collapse,  we  are  not  concerned  about  the  details  of  the  SSBs  geometry.  Also  for  larger  domains  and 
coarser  resolution  we  can  simulate  the  SSB  using  a  subgrid  scale  (SGS)  model.  The  SGS  model 
includes  the  SSB  as  a  point  source  of  momentum,  e.g.  Sirviente  and  Patel  (1997) .  By  using  the  SGS 
model  as  a  lateral  BC  we  can  extend  the  simulation  results  to  regions  remote  from  the  SSB. 

WORK  COMPLETED 

Scaling  analysis  and  an  examination  of  the  limitations  on  computational  domain  size  led  us  to 
parameterize  the  submarine  as  a  momentum  source  (i.e.,  a  jet)  at  the  simulation  boundary.  Attempting 
to  simulate  the  submarine  as  a  body  of  finite  size  would  have  restricted  us  to  a  domain  too  close  to  the 
submarine.  We  are  interested  in  the  mid-range  to  far  field  jet  of  the  submarine  and  its  interactions  with 
the  ambient  internal  wave  field. 

We  first  simulated  a  jet  in  an  unstratified,  homogeneous  fluid.  The  results  were  compared  with 
analytically  derived  similarity  profiles  and  laboratory  data.  We  then  simulated  a  jet  in  linearly  stratified 
fluid.  The  ratio  of  vertical  to  horizontal  kinetic  energy  was  calculated  for  several  stratifications.  The 
results  were  similar  to  those  found  in  numerous  laboratory  experiments  and  other  numerical 
experiments. 

Finally,  we  have  started  inserting  a  jet  into  a  realistic  simulation  of  the  ambient  temperature,  salinity 
and  velocity  in  the  vicinity  of  the  Alabama  Alps  in  the  Gulf  of  Mexico  (GoM).  This  simulation 
includes  full  Alabama  Alps  bathymetry  and  is  initialized  with  observed  temperature  and  salinity 
profiles.  The  Alabama  Alps  significantly  affect  the  flow  in  their  vicinity.  The  flow  over  the  Alps 
generates  hydraulic  jumps,  trapped  lee  waves,  large  amplitude  internal  waves  and  vortices  that 
propagate  downstream.  Also,  temperature  and  salinity  are  mixed  and  turbulent  boundary  layers  form 
behind  the  alps  for  the  higher  velocity  flows. 
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RESULTS 


HOMOGENEOUS  FLUID 

The  first  experiments  involved  a  round  axisymmetric  jet  flowing  into  a  homogeneous  fluid.  The 
parameters  for  this  experiment  are  given  in  Table  1. 


Table  1:  Parameters  used  in  homogeneous  jet  simulations. 


Lx(  m) 

Lv(m) 

Lf  m) 

zk(m) 

Zly(m) 

zfz(in) 

Uo  (m/s) 

D(  m) 

At  (s) 

480 

160 

160 

1 

1 

1 

0.2 

10 

1 

Using  Uo  and  D  as  velocity  and  length  scales,  the  Reynolds  number  is  defined  as: 

U,D 


ReD  — 


v 


(1) 


In  this  study  two  Reynolds  number  simulations  are  performed,  Reo=400  and  Re«  =2000.  Figure  1 
contains  horizontal  and  vertical  cross-sections  of  the  flow  evolution  for  each  Rqd  at  t= 10800s.  The  Re^ 
=400  simulation  appears  to  be  laminar  with  a  nicely  formed  vortex  pair  at  the  jet  leading  edge.  In 
contrast,  the  Re^  =2000  simulation  appears  to  be  turbulent,  spreading  faster  and  traveling  slower 
downstream  than  the  laminar  jet.  The  turbulent  jet  has  larger  velocity  gradients  than  laminar  flow, 
resulting  in  a  higher  rate  of  energy  dissipation  and  ambient  fluid  entrainment. 


Re  =  400 


Re  =  2000 


Horizontal 
cross  section 


Vertical 
cross  section 


Figure  1.  (Upper)  Horizontal  and  (lower)  vertical  cross-section  of  U  for  each  Re.  Both  plots  are  at  t 
=  10800s.  The  Reo=400 jet  is  laminar  with  a  vortex  pair  formed  at  the  leading  edge  of  the  jet.  The 
Reo=2000 jet  is  turbulent  with  large  velocity  gradients  seen.  Due  to  large  velocity  gradients  the 
turbulent  jet  dissipates  more  energy  than  the  laminar  jet  and  does  not  extend  into  the  domain  as 

quickly  as  the  laminar  jet. 


As  each  jet  progresses  downstream  the  initial  top-hat  jet  becomes  Gaussian  in  shape.  Analytic  results 
(e.g.,  Pope  2000,  Tennekes  and  Lumley  1972)  demonstrate  that  the  jet  will  become  Gaussian  and  will 
be  self-similar  when  scaled  by  the  centerline  velocity,  Uc,  and  jet  half  width,  r/2.  The  scaled  velocity 
profiles  are  shown  in  Figure  2,  and  fit  well  to  a  Gaussian  curve  after  x/D=\0.  The  profiles  are  the  result 
of  averaging  over  50  snapshots  of  the  model.  The  turbulent  profiles  do  not  fit  the  Gaussian  as  well  as 
the  laminar  jet,  possibly  because  the  averaging  period  was  too  short  to  minimize  turbulent  fluctuations. 
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Figure  2.  Normalized  horizontal  (xy)  plane  profiles  of  U  for  Re  [,=400  (left)  and  Re  i, =2000  (right).  A 
top  hat  shape  is  seen  at  the  jet  origin  x/D=0.  The  profiles  become  self  similar  (solid  black  line)  by 
x/D=5  when  scaled  by  centerline  velocity,  U„  and  jet  half  width,  r/2.  Noise  in  the  Re  [,=2000 jet  may 
be  a  result  of  turbulent  fluctuations  that  have  not  been  removed  via  ensemble  averaging. 

LINEARLY  STRATIFIED  FLUID 

In  the  second  set  of  simulations  the  jet  flows  into  a  linearly  stratified  fluid.  In  stratified  flows  buoyancy 
forces  exist  that  can  affect  the  flow.  The  Froude  number  is  now  defined,  which  is  the  ratio  of  inertial  to 
buoyancy  forces: 


(2) 


stratification.  Stratification  inhibits  the  vertical  spreading  of  the  jet  due  to  the  buoyancy  forces  (Figure 
3),  and  the  jet  becomes  asymmetric  compared  to  the  jet  in  homogeneous  fluid.  For  F=1  the  outer  jet 
edge  is  wavy  and  demonstrates  alternating  positive  and  negative  x-direction  flow  suggesting  a  weak 
von  Karman  vortex  street  has  formed. 

The  vertical  extent  of  the  jet  does  not  increase  with  decreasing  Froude  number  as  anticipated.  This 
might  be  explained  by  the  relative  short  time  of  the  simulation  whereby  any  noticeable  vertical 
difference  have  not  formed  yet  or  it  may  be  an  artifact  of  the  contouring.  The  impact  of  stratification 

can  be  seen,  however,  in  plots  of  the  volume  averaged  horizontal,  (Eh),  and  vertical,  (Ev) ,  kinetic 

energy  (Figure  4).  As  expected,  the  vertical  kinetic  energy  decreases  and  the  horizontal  kinetic  energy 
increases  with  decreasing  Froude  number,  increasing  stratification.  This  result  is  in  agreement  with 
analytical  studies  (e.g.,  Lilly  1983)  and  numerical  experiments  (e.g.,  Riley  and  de  Bruyn  Kops  2003). 
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Horizontal  cross  section 


Vertical  cross  section 


Figure  3.  (Left)  Horizontal  and  (Right)  Vertical  slices  ofU  for  F=l,2,4,  Reo=2000,  t=l 0800s.  The 
horizontal  cross  section  shows  a  wavy  pattern  with  oscillating  velocity  for  the  strongly  stratified  1=1 
case,  reminiscent  of  a  vortex  street.  There  is  not  noticeable  different  in  the  vertical  heights  of  each 
jet.  This  is  likely  a  result  of  the  jet  not  having  enough  time  to  fully  develop. 


Figure  4.  Evolution  of  (Eh)  (left)  and  {Ev)  (right)  for  each  Froude  number  and Re=2000.  The 

homogeneous  flow  is  plotted  as  a  black  line  for  comparison.  A  clear  trend  is  seen  where  as  the 
stratification  becomes  stronger  (i.e.,  F  decreases),  the  jet  decays  slower  and  has  reduced  vertical 
kinetic  energy.  This  indicates  the  flow  becomes  less  3D  (or  more  2D)  with  increased  stratification 
and  vertical  decoupling  of  the  flow,  in  agreement  with  theoretical  (e.g.,  Lilly  1983)  and  numerical 

studies  (e.g.,  Riley  and  de  Bruyn  Kops,  2003). 
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REALISTIC  STRATIFICATION/SIMULATIONS  IN  THE  GULF  OF  MEXICO 


The  motivation  for  the  work  in  this  section  is  to  simulate  the  wake  of  a  submarine  near  variable 
topography  and  in  realistically  stratified  water.  The  Alabama  Alps,  which  is  a  submerged  ridge  in  the 
GoM,  were  chosen.  Common  to  all  simulations  in  this  section  are  the  initial  temperature,  salinity, 
density  and  buoyancy  frequency  profiles  shown  in  Figure  5. 


TlaC) 


5  (psu] 


Figure  5.  Initial  temperature  (T),  salinity  (S),  density  ( <Jd),  and  buoyancy  frequency  (N).  These 
profiles  are  a  smoothed  average  of  observations  taken  at  12  stations  around  the  Alabama  Alps. 

A  series  of  simulations  of  the  jet  were  conducted  to  create  initial  conditions  to  insert  into  the  full  GoM 
domain.  The  parameters  for  these  simulations  are  given  in  Table  2.  This  resulting  initial  condition  will 
then  be  inserted  into  the  results  of  a  GoM  simulation  that  was  run  to  a  fully  developed  state  without  a 
jet.  All  simulations  were  conducted  with  the  centerline  of  the  jet  at  mid-depth  (50m). 


Table  2:  Parameters  for  GoM  initialization  simulations. 


Lx(  m) 

Lv(  m) 

Lz  (in) 

Nx 

Ny 

Nz 

dx(m) 

M  m) 

Az{ m) 

At(  s) 

1024 

2048 

100 

256 

512 

101 

4 

4 

1 

1 

The  intent  of  the  simulations  is  to  study  the  interaction  of  the  intermediate  and  far  field  of  the  jet  with 
an  ambient  IW  field,  particularly  large  amplitude,  nonlinear  internal  waves.  This  prohibits  the  detailed, 
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near  field  modeling  of  the  SSB  because  at  it  would  require  very  high  resolution  in  the  vicinity  of  the 
SSB.  Even  with  grid  refinement  techniques  it  is  not  computationally  feasible  to  simulate  both  the  near 
and  the  far  field.  Thus,  the  simulation  of  the  SSB  starts  downstream  from  the  body  itself. 

It  is  also  desirable  to  start  the  simulation  before  the  wake  begins  to  collapse.  The  onset  of  wake 
collapse  begins  at  time  tc=2n/3N  (Kaulum,1972).  At  z=50m  (where  the  jet  is  centered),  A=0.0148 
rad/s,  resulting  in  tc  =141  s. 

For  these  simulations  the  speed  of  the  SSB,  V,  is  assumed  to  be  lOm/s  (~20knots).  Thus,  the 
downstream  position  at  the  onset  of  wake  collapse  is:  xc=Vtc  =  1410m.  The  downstream  starting 
position  of  the  SSB  is  arbitrarily  set  at  x=500m,  which  puts  the  jet  well  before  the  start  of  the  collapse 
stage,  and  unstratified  relationships  can  be  used  to  determine  the  jet  maximum  velocity  and  diameter. 

For  unstratified  conditions,  Kaulum  (1972)  provides  the  following  relationship  to  obtain  the  wake 
diameter,  D,  at  doenstream  position,  x, 


D  f  x  Y 

—  =  1.251  — I 

D0  { D0) 


(3) 


where  D0  is  the  initial  wake  diameter  and  n= 0.22  is  an  empirically  obtained  exponent.  Substituting 
x=500m  into  (3)  yields  a  jet  diameter  D  =  29.6m,  which  is  rounded  to  30m  for  these  simulations. 

A  time  series  of  horizontal  cross  sections  of  the  downstream  velocity  for  the  jet  simulation  in  GoM 
stratification  is  shown  in  Figure  6.  By  hour  2  the  jet  has  penetrated  across  half  the  domain  and  strongly 
impacts  the  middle  third  of  the  domain.  After  eight  hours  the  jet  extends  across  the  domain.  The  jet 
takes  up  a  significant  portion  of  the  vertical  extent  (Figure  7),  as  it  is  30m  in  diameter  vs.  101m 
domain  depth.  The  jet  magnitude  decays  as  it  continues  to  propagate  across  the  domain  by  hour  10. 


t=0hr  t=2hr  t=4hr  t=6hr  t=8hr  t=l  Ohr 


Figure  6.  Time  series  of  horizontal  cross-sections  of  V  through  jet  center  (z=50m).  Color  scaling  is 
the  same  as  in  Figure  7.  The  jet  progresses  into  the  domain  and  becomes  wavy  by  t=10hr, 
reminiscent  of  a  vortex  street  from  the  wake  of  an  object. 

A  horizontal  temperature  cross  section  (Figure  8)  shows  the  horizontal  variability  induced  by  the  jet. 
Waves  are  seen  on  the  leading  edge  and  trailing  off  the  front  of  the  jet  at  an  angle.  Vortices  begin  to 
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form  at  hour  7  and  by  hour  10  a  von  Karman  vortex  street  has  formed.  Internal  waves  are  forming  in 
the  vertical  temperature  cross  section  (Figure  9).  The  waves  appear  to  be  significatntly  nonlinear  and 
the  greatest  wave  formation  is  in  the  center  of  the  domain  under  the  head  of  the  jet.  An  isosurface  of 
temperature  in  the  horizontal  cross  section  through  the  center  of  the  jet  (Figure  10)  clearly  shows  the 
creation  of  nonlinear  internal  waves  in  front  of  and  around  the  jet  as  well  as  propagation  of  the 
nonlinear  internal  waves  away  from  the  jet. 


-0.1  0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 


Figure  7.  Time  series  of  vertical  cross  section  of  downstream  velocity,  V,  through  the  jet  center. 
Units  on  color  bar  are  in  m/s.  The  patchiness  in  the  velocity  downstream  is  caused  by  the  oscillation 
of  the  jet  at  it  progresses  through  the  simulation  domain  (Figure  6).  Note  the  jet  takes  up 
considerable  space  as  it  is  30m  in  diameter,  whereas  the  vertical  extend  of  the  domain  is  101m. 
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Figure  8.  Temperature  contours  through  the  center  of  the  jet.  A  leading  wave  can  be  seen  in  the 
temperature  field  at  t=2hrs.  Also  there  seems  to  be  a  wave  generated  by  the  jet  that  propagates  at  a 
near  45  degree  angle  away  from  the  jet.  By  t=10  swirls  of  temperature  indicating  that  a  von  Karman 

vortex  street  has  formed. 


Figure  9.  Temperature  contours  at  t=10hrs.  The  structure  seen  in  the  fluctuations  near  the  center  of 
the  domain  are  nonlinear  internal  waves  formed  by  the  jet. 
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Figure  10.  Isosurface  of  temperature  at  10 hr s  on  a  horizontal  cross  section  through  the  centerline 
of  the  jet.  The  jet  creates  vertical  motions  that  displace  temperature.  Wave  motions  can  also  be  seen 

propagating  away  from  the  jet. 

Simulations  of  the  flow  around  the  Alabama  Alps  have  been  conducted  with  the  temperature  and 
salinity  profiles  discussed  above  as  initial  conditions  and  a  barotropic  eastward  current  as  forcing. 
Experiments  were  conducted  with  three  forcing  speeds,  5  cm/s,  10  cm/s  and  50  cm/s.  In  the  5cm/s  case 
the  flow  is  blocked  by  the  Alps.  The  flow  is  hydraulically  controlled  in  the  10  cm/s  case  and  in  the  50 
cm/s  experiment  the  flow  is  supercritical,  with  respect  to  the  Froude  number,  and  nonlinear  internal 
waves  form  behind  the  Alps  and  disperse  downstream  (Figure  11).  This  work  has  been  submitted  to 
Geophysical  Research  Letters  (Gallacher,  et  al.,  2008). 
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Figure  11:  Streamlines  for  (a)  the  5cm/s  case;  (b)  the  lOcm/s  case  and  (c)  50cm/s  case.  Streamlines 
bend  around  the  alps  for  the  5  and  lOcm/s  simulations,  indicating  flow  blocking,  whereas  the  fluid 

flows  over  the  alps  in  the  50cm/s  simulation. 

Although  the  forcing  velocity  is  consant  in  time  and  space,  the  flow  resulting  from  the  interaction  with 
the  Alps  is  quite  complex  with  significant  spatial  and  temporal  variability.  This  can  be  seen  in  the 
kinetic  energy  dissipation  rate  shown  in  Figure  12.  Also,  the  difference  in  the  kinetic  energy 
dissipation  rate  for  the  10  cm/s  case  versus  the  50  cm/s  case  is  substantial. 


Figure  12.  Kinetic  energy  dissipation  rate  for  U=10cm/s  (left)  and  50cm/s  (right).  The 
50cm/s  simulation  is  more  turbulent  and  generates  a  larger  wake  with  more  mixing. 


IMPACT/APPLICATIONS 

This  work  will  help  to  determine  the  importance  of  and  the  requirements  for  nonhydrostatic  forecast 
systems  for  naval  applications.  The  scales  and  features  which  will  require  nonhydrostatic  simulation 
are  being  assessed. 

RELATED  PROJECTS 

The  NRL  project  Autonomous  Characterization  of  Environmentally  Induced  Non- Acoustic  Noise  and 
the  Adaptation  of  Multi- Sensor  USW  Networks.  (6.2,  Undersea  Warfare)  is  related  to  this  project 
because  it  involves  nonhydrostatic  modeling  of  the  SW06  experimental  area  and  time  and  comparison 
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with  measurements  taken  during  SW06.  Components  of  SW06  are  funded  through  this  ONR  NLIWI 
DRI. 
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